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NOTATION 


Dotted  variables  have  dimensions;  undotted  variables  are 
dimensionless  (reference  length  is  U,s/g).  Symbol  ~  stands  for 
Fourier  transform.  Superscripts  s,  b,  l  stand  for  free-surface . 
body  and  waterline  effects,  respectively. 


b1 

B' 


s  b  „  t 
Ci  ,  Ca  >  C a  ,  C 


n 

a  >•  a 


f(x,z) 

Fn 

g 


G 

_ 

GrjG 

h(x) 


Ii  ,  la  ,  la  >  la 


L' 


L(a,v) 

m(z) 


half  beam 
beam 

first  and  second  order  cosine  spectrum 
functions 

hull  shape  function 

length  Froude  number 

acceleration  of  gravity 

Green  function 

parts  of  the  Green  function 

waterline  shape  function 

auxiliary  function 

auxiliary  functions 

half  length  of  the  ship  (<<  »  ^'g/lP8) 

length  of  the  ship 

auxiliary  function 

cross-section  shape  function 
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s  b  t 
M  ,M.,M 
J 

s  b  t 
P  2 5  P  2 i Ps 


s  b  t  s  D  l 
q  %q  oQ  ,q  ,q.  ,q 

<3  P  r1,P  P 


R 


s  b  l 
Si  jSa^SajSg^Sa 


t' 

U' 

x' ,y‘ ,z« 


auxiliary  functions 

auxiliary  functions  equivalent  to 
pressures 

auxiliary  functions 
wave  drag 

submersed  area  in  the  center  plane 

first  and  second  order  sine  spectrum 
functions 

draft 

ship  velocity 

cartesian  coordinates  attached  to  the 
ship 


a 

P 

ft(z) 

C  «  b'/t* 

€  =  b ' g/U ' 9 

>®a  »$a 

w  ,sw  bw  lw  w 
><Pa  ,«Pa  ,«P*  ,«P 


longitudinal  and  transversal  wave  numbers 

root  of  aa  -7-0 

Dirac  function 

slenderness  parameter 

slenderness  parameter 

first  and  second  order  potentials 

free  wave  potentials 
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V 

r 


(as  +  ?)K 

[(a-v)a  +  (i?-T)ap 


v 

P  i 

p  =  (v2  +  Ta)^ 

T 

Ci,cLc*,ctc> 


.  sw  rbw  Aw  rw 

»  X  j  »  3  3  t*a  ,  U 


auxiliary  functions 

coefficient  of  artificial  viscosity 

integration  variable 

root  of  va  -  p  =  0 

auxiliary  function 

integration  variable 

first  and  second  order  free  surface 
elevation 

first  and  second  order  free  waves  profile 
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ABS TRACT 

The  present  work  investigates  second  order  terms  of  the 
expansion  in  the  beam/length  ratio  of  the  potential  of  flow  past 
a  thin  ship.  The  ship  is  in  a  steady  translational  motion,  in 
water  of  infinite  depth  and  infinite  lateral  extent,  within  its 
central  plane  of  symmetry.  Heave  and  trim  are  not  taken  into 
account . 

The  expressions  suggested  by  Maruo  (1966)  for  the  second 
order  potentials  are  taken  as  the  starting  point.  The  second 
order  potential  is  split  into  the  free-surface,  body  and  water¬ 
line  integral  corrections.  The  associated  sine  and  cosine 
spectrum  functions  are  derived  systematically  by  using  double 
Fourier  transforms.  Computation  of  the  value  of  a  second  order 
spectral  function  for  a  given  transversal  wave  number  is  shown 
to  require  six  integrations  in  general  and  two  integrations  in 
the  case  of  a  simple  analytical  hull  shape  function.  Simplifica¬ 
tions  of  spectral  functions  in  the  case  of  separable  hull  func¬ 
tions  and  symmetrical  ships  are  discussed. 

The  second  order  spectral  functions  are  computed  in  detail 

for  the  particular  case  of  a  ship  with  parabolical  waterline  and 

rectangular  frames.  Preliminary  numerical  results  for  the  second 

order  sine  spectral  furv  Uons  for  the  draft/lgngth  ratio  equal  to 

0.15  and  for  the  Proude  numbers  P  *■  0.316  and  P  &  0.200  are 

r-  n 

presented  and  discussed. 
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l.  INTRODUCTION 


We  consider  the  generation  of  waves  and  the  associated  wave 
resistance  of  a  thir.  ship  moving  at  constant  speed  in  still  water 
of  infinite  depth  and  infinite  lateral  extent.  The  usual  linear¬ 
ised  theory  (Michell's  theory)  is  known  to  be  a  first  order  solu¬ 
tion  in  the  beam/length  ratio  (Peters  and  Stoker,  1967).  The 
general  expressions  of  the  second  order  potential,  and  the  associa 
ted  wave  resistance,  have  been  derived  in  a  complete  form  by 
Wehausen  (196?)  and  Maruo  (1966).  The  aim  of  the  present  work 
is  to  derive  a  method  of  computing  the  second  order  effects, 
which  has  become  possible  with  the  availability  of  fast  electronic 
computers. 


Duch  computations  are  extremely  valuable  for  the  further 
advancement  of  the  theory  of  wave  resistance  and  its  application. 
They  may  result  in  a  more  accurate  -  valuation  *jf  wav?  r**Utane* 
of  ships  tm  those  case;  in  which  the  linearised  .  theory  has  been 
found  to  be  unsatisfactory .  steend  computations 


may  also  help  ir.  dec  id  inf  whether  a^ac  of  the 
foctb  tan  be  n^fi- ,_„ed  in  sruer  to  arrive  at  i 


gc-cend  order  «f« 
neons i. stent ,  but 


simpler,  onsviatior.s.  *\trtht  rarer#*  by  tbssparlnf  measured  with 
cbmpute^  wave  sport ra  it  Is  possible  to  find  out  whether  second 
order  effect  a  or  virebus  ef  fecte  ere  re-ep.bneible  f#r  the  dir*, 
crepanty  between  the  first  order  solution  and  experiments. 
Ktn-ally*  derivation  of  shapes  «f  little  wave  resistance  can  be 


improved  by  accounting  for  nonlinear  effects. 
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The  second  order  wave  resistance  has  been  computed  by 
Eggers  (1970)  for  a  parabolical  ship.  The  present  work  parallels 
somehow  Eggers  fundamental  studies,  but  differs  in  method:  the 
second  order  free  waves  potential  is  derived  by  straightforward 
double  Fourier  transform  and  the  wave  resistance  coefficient  is 
computed  from  the  free  waves  spectra. 

The  influence  of  trim  ancl  heave  is  not  taken  into  account 
here,  but  the  method  can  be  used  in  order  to  incorporate  these 
effects  too. 


II.  GENERAL  FORMULAE  FOR  FIRST  AND  SECOND 

ORDER  POTENTIALS  AND  FOR  WAVE  RESISTANCE 

We  consider  a  symmetrical  ship  of  length  L'  =  2^',  beam 
B'  =  2b'  and  submersed  area  (in  the  center  plane,  beneath  the 
unperturbed  level)  S',  moving  at  constant  speed  U'  in  the  posi¬ 
tive  x'  direction  (Fig.  1)  within  its  plane  of  symmetry.  The 
axes  x'  and  y!  lie  in  the  horizontal  plane  of  the  unperturbed 
free-surface  while  z?  is  vertical  and  positive  upwards,  z*  = 
-t'(x')  describes  the  contour  of  S'. 

The  variables  are  first  made  dimensionless  with  respect  to 
U'  and  U'2/g,  i.e. 

x,y,z  =  x ' g/U ' 3 ,  y ' g/U ' 8  ,  z'g/TJ'2  ;  cp  =  cp'g/U'3 

[2.1] 

C  =  C'g/U'a  ;  t  =  t'g/U'8  ;  l  =  t'g/U'8  ;  e  =  b'g/U'3 
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where  cp '  (x1  ,y '  }z ' )  is  the  velocity  potential  (the  velocity  is 
its  positive  gradient)  and  C 1 ( x * , y ^ )  is  the  free  surface  elevation. 

The  equation  of  the  hull  surface  is  as  follows 


y  =  *  €f(x*z)  ( U!  <  is  z  >  -t)  [2.2] 


We  consider  the  translational  steady  motion  of  the  ship  which  is 
kept  otherwise  fixed  with  .aspect  to  the  coordinate  system. 

Assuming  that  e  =  o(l)  while  l  =  0(1)  and  t  =  0(1)  the 
potential  is  expanded  in  a  perturbation  series  as  follows 


*(*,y,ai6,*,...)  ~  -x  +  ecpx(x,y,z;^,...)  +  escp3  (x,y,  zj-t, . . . )  +  .. 


[2.3] 

For  the  -  ake  of  completeness  we  give  here  the  expressions  of  cpx 
and  c p3  following  Maruo  (1966)  formulation: 


«Pi(x,y,z)  ff  ^G(x,y,zj  x,0,z)  dxdz  [2.4] 


s  b  K> 
cpa  =  cp8  +  cp3  +  cp3 


[2.3] 


¥8(x,y,z) 


4m 


Pa(x,y)  G(x,y,z;  x,y,0)  dxdy[2.6] 


-00 
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cpa(x,y,z)  =  ~ JJp a(x,z)  G(x,y,x;  x,0,z)  dxdz  [2,7] 


-b  f  p» 

-t 


Pa (x)  G(x,y,z;  x,0,0)  dx  [2.8] 


s  b  ^ 

where  cp3,  coa  and  co2  are  the  free-surface,  the  body  and  the  water¬ 
line  integral  second-order  corrections,  respectively,  and  G  is  the 
Green  function.  The  pa  functions  have  the  following  expressions: 


S/~ 


dcPx  dac0;L  dcp-L  dacoi  dcpi  dacpx 


pa(x,y)=  3-^  —  +  2 


dx* 


3y  9xdy  dx 


dcp1  Sacpx  dcpi  da  cpx 


•'  2 


dz  dxdz  dx 

dxddz 


(x=x,  y=y,  z=0)  [2.9] 


1,, — ,  a  ,  St<\ 

Pa(x,z)=  -  jj  (f-Jj) 


dcpi 

> 

dz' 


( x=x,  y=0,  z=z)  [2.10] 


e°(x>=  5J  U 


(x=x,  y=0,  z=0)  [2.11] 


The  Green  function  is  given  by 

1 


G(x,y,z;  x,y,z)  = 


[(x-x)a  +  (y-y)3  +  (z-z)3]2 


[(x-x)a  +  (y-y) 3  +  (z+z)ar 


+ 


[2.12] 
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where  Gr  is  regular  in  the  lower  half  plane  2  <  0  and  G  satisfies 
the  linearized  free-surface  and  radiation  conditions. 


We  shall  use  in  the  sequel  the  following  Fourier  transforms 


CO  00 

*  -  5F  /  ax  / 


dy  *(x, y,z)  [2. 13] 


go  a 

I  aaf 


dp  $(a,0,z)  e 


-i(ax+J3y) 


[2.14] 


where  a  and  /3  are  the  longitudinal  and  transversal  wave  numbers, 
respectively. 

The  Fourier  transform  (FT)  of  G  has  the  following  expression 
(see,  for  instance,  Timman  and  Vossers,  1955) 


~  hr 

G  =  (—  ( 

r  L 


e~Y | z-z |  _  eY(z+z) 


sy(z+z) 


a  -  y 


-  i|ja 


i(ax+/3y) 


[2.15] 


where  y  =  (a®  +  jS8)2  is  real  and  positive  for  a,  p  reals  and  u 
is  a  positive  coefficient  of  artificial  viscosity  which  is  let  to 
tend  to  zero  (u  is  introduced  for  an  easy  account  of  integration 
paths  in  the  complex  plane).  The  three  terms  of  G  [2.15J  are  the 
FT  of  the  corresponding  terms  of  G  [2.12], 
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Tne  FT  of  the  different  potentials  are  obtained  by  replacing 
G  by  G  in  [2.1;]  -  [2.8].  We  are  interested  here  mainly  in  the 
potentials  of  the  free  waves  for  x  ~  -®.  Only  part  of  G,  namely 


£W_  2  ei(ax  +  £y) 

as  -  y  -  iua 


[2.16J 


contributes  to  these  potentials.  We  can  write,  therefore. 


5?  -  57-  f  f  ~  GW(x,0)zj  z)  dxdz  (x  -  -")  1.2.17] 

JSJ 

...  .  ~sw  ~bw  tw 

and  similar  expressions  for  cpa  ,  cp8  and  cp8  . 

To  simplify  notations  we  define  the  following  functions 


Ii(o.,/3)  =  JJ*  ;  eYZ  e‘aX  dxdz  [2.18] 

S  9x 

Ia(a,j3)  =  |  J  dx  J’  dy  pf(x,y)  [2.19] 

-08  _00 

la ''a,  A)  =  J’J*  P9  (x,  z)  e7Z  eiaX  dxdz  [2.20] 

S 

l 

la  (a,  A)  =  J’  pa(x)  eiaX  dx 

-l 


[2.21] 
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The  free  wave  potential  cpi  can  be  obtained  from  cpx  [2.17J 
by  inversion  [2.14].  The  free  waves  potential  results  from  the 
contribution  of  the  residues  at  the  poles 

a  =  ±  ap  ;  ap  =  j^l  +  (1+4j33)^  J  [2.22J 

of  the' quotient  l/(a  -y)  of  G  1.2.16  j  for  integration  over  a. 

The  presence  of  u  in  [2.16]  requires  circumveting  these  poles 

from  below.  The  residue  of  the  quotient  is  2-rria  /(2a  8-l)  while 

P  v  P  ’ 

■v(ap,j3)  =  y(-ap,£)  =  ap2. 

Hence,  the  free  waves  potential  cp^  may  be  expressed  with 
the  aid  of  [2.16],  [2.17]  and  [2.18]  as  follows 


00 

<P?(x,y,z)  =  ^-i^(ap,0)  -  Ii(-ap,/3) 


cos  a  x  - 
P 


Ii  (ap,P).  +  Ii  (-ap,0) 


\  a  a 
a  x)  - e  1 

p  J  2a  3 -1 


sin  a  x  >  — — —  e  p  cos  jSyd/3 

[2. 23] 

sw  bw  tw  r 

the  expressions  of  cp8  ,  cp8  and  cp8  being  similar  to  [2.23J  with 

Ii  replaced  by  I8 ,  I8  and  I8  [2.19]  -  [2.2l],  respectively.  Cos  /By 
has  been  substituted  for  in  [2.23]  because  of  the  sym¬ 
metry  of  the  flow  with  respect  to  y  *  0.  Equations  [2.18J  -  12.21J 
show  indeed  that  lj(a,0)  -  I ^  ( -a , /S)  is  purely  imaginary  while 
Ij(a,^)  +  Ij(-a,/3)  is  real  (j  »  1,2),  so  that  cp^  are  real. 
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The  equation  of  the  free -surface  has  an  expansion  similar 
to  L2.3] 


C(x,y;e,t,...)  ~  eCj (x,y; l, . . . )  +  e8Ca (x,y;t, . . . )  +  ... 

1.2.24] 

The  equation  of  the  free  waves  profile  is  obtained  from  the 
potential  as  follows 


w 


C?(x,y)  = 


(x,y,0) 

5x 


[2.25'] 


.w 


,  sw  ,bw  M  a^w(x,y,0) 


C."(x,y)  =  Ci  +  +  £.  = 


ax 

*PaW(x,y,o) 

— 


ax 


[2.26] 


w 


It  Is  worthwhile  to  mention  here  that  £,  differs  from  q  by  local  terms 
which  decay  algebraically  like  (x3+y8)*'a  for  x  -  while  the 
local  terms  in  Ca  are  of  order  (x8+y8)”*. 

It  is  seen  from  [2.25]  and  [2.23J  that  one  can  write  for 
the  free  waves  profile 


C?(x,y) 


00 

-  r  /  t* 


a 

(0)  sin  a  x  i  CiO)  cos  a  x]  008 


t2.27] 
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where  the  sine  and  cosine  spectrum  functions  Si ( 0)  and  Ci(£) 
are  given  by 


Si(jB)  -  iCi  O)  =  2ili(a  ,  [2.28] 

r 

■-v..-  bw  tw 

Similar  expressions  are  obtained  for  Ca  >  Ca  and  Ca  with  S*  3 

Ci  replaced  by  Sa,  Ca,  Sa,  Ca,  Sa,  Ob’,  respectively,  in  [2.27j 
and  Ji.  by  la,  la  and  la  in  [2.28], 

Finally,  the  dimensionless  wave  drag  R  =  R'ga/pU'4  can  be 
written,  on  the  basis  of  the  first  and  second  order  potentials, 
as  follows 

•:  .  0O 

\  r  a  8 

R-~  /  [(eSi  +  e8Sa)8  +  (eCv  +  e8Ca)8]  - d/3  [2.29] 

2lr  -»  2a  8 -1 

P 

Expending  R  in  an  e  series  yields 


1?:.  R  ~  e#Ri  +  e3Ra  +  0(e4) 

Ji-  ■  > 

where 


a 

+  C18)  - E —  d/3 

2a  8-l 
P 


Ra  =  i  J  (SiS.  +  Ci  Ca  )  - dj3 

*  -•  2a  8-l 

P 


[2.30J 


[2.31] 


[2.32] 
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R  [2.29]  is  always  positive,  whether  R  [2.30]  may  become 
negative  if  the  second  order  effects  are  exceedingly  large. 

Using  [2.29]  is,  however,  inconsistent  because  it  includes 
contributions  0(e4). 

All  the  formulae  of  this  section  are  well  documented  in 
the  literature  (for  instance  in  Wehausen,  1973)  and  have  been 
given  here  for  the  sake  of  completeness  and  clarity  of  notation. 

The  purpose  of  the  present  work  is  to  derive  the  second 
order  spectrum  functions  Ss  ( /3)  and  Ca(|S)  for  ships  of  given 
shape. 


III.  THE  COMPUTATION  OF  THE  SECOND  ORDER 
FREE  WAVES  SPECTRUM  FUNCTIONS 

Equations  [2.191  -  [2 . 21 ]  render  the  functions  needed  for 
the  computation  of  the  spectral  functions  in  terms  of  the  first 
order  solution  cpj. .  Considerable  simplifications  in  number  of 
integrations  and  smoothness  of  integrands  may  be  achieved  how¬ 
ever,  by  further  elaboration  of  the  integrals  of  [2.19]  -  [2.21]. 

1.  The  Complete  Expression  of 

Equations  [2.16]  and  [2.17]  render  the  expressions  of  the 
FT  of  the  first  order  free  waves  potential.  Hie  second  order 
terms,  however,  are  based  on  the  complete  expression  of  ®x,  i.e., 
on  the  local  as  well  as  the  far  waves  term.  It  is  useful,  there¬ 
fore,  to  write  ^  and  its  derivatives  in  an  explicit  form.  From 
[2.4]  and  [2.15]  we  have 
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-y|z-z[  7(2+2 ) 


-  e 


1- 


0  y ( z+z )  ,  . 

2 e  K  'l  iax 

-  )  e  dxd2 

as  -7-iUa 


[3.1] 


We  shall  consider  here  smooth  ship  shapes,  the  highest 
conceivable  singularity  of  f  being  associated  with  a  flat  bottom, 
where  f  has  a  finite  jump  for  z  =  -t.  It  is  seen  from  [3.l]  ^that 
cpx  and  (Bcpj/dz)  are  then  continuous  for  z  =  -t  and  (d*cpi/dza)  has 
a  finite  discontinuity  at  z  =  -t.  This  observation  will  be  used 
in  the  sequel. 

cpa  [2.6]  and  cpj  [2.8]  are  based  on  the  equations  of 
cpi(x,y.O)  and  its  derivatives.  For  z  «  0  we  obtain  from  [3»l] 
the  following  relationships 

5i(a,0,o)  =  0)  -  -  t3.?J 


T3T  *  '  ia{pl  *  “57  "  -  l&''> 


a  <Pi 


(z  -  0)  [3.3] 


and  similarly  for  higher  order  derivatives. 
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cpa  t2.7j  is  based  on  dcpi/dx  and  dcpx/dz  along  S  (y=0)„  For 


the  FT  of  these  components  we  obtain  from  [3.l] 


dePi 

~5x 


=  -  iacpi 


(y  =  o) 


[3.4] 


dcpv 

"TjT 


7  =  ^// |  [-  [e-'^'sUn  (z-I)  +  e^->]  - 


2;ey<z+*'>  1  lax  -- 

— -  >  e  dxdz 

a*-v-ltia 


[3.5] 


2.  The  Computation  of  S>  and  C»  . 

We  consider  first  the  free  surface  second  order  effect.  The 
function  l!(a,p)  [2.19]  is  proportional  to  the  FT  of  ps(x,y) 

[2.9].  Using  the  convolution  theorem  in  itB  following  symmetrical 
form 


m  m 

m(x^yM(xIy)  •  ^/dvf  dT  C«(a-v,i3-T)?(v,T)+w(v,t)T(a-v,^-T)] 

13-63 


where  •  and  9  are  arbitrary  functions  having  FT  and  v,t  are 
integration  variables,  we  may  write  with  the  aid  of  [2.9]  and 
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00  00 

!*(<*,£)  =  j  f  dv /  dT  hS(aJjSJs;,T)^(v,T,o)^(a-vJ^-T,0)  [3.7] 

_C0  .00 


with 

hs(a,/3,v,t)  =  3av(a-v)+2aT  (0-t  )+v(  /3-t  )*+(a-v)“v8-2aMa (a-v)*- 


-  v(a-v) 

[v3+(a -v)a  ] 

[3.8] 

By  substituting  in  [3.7]  the  expressions  of  [3.2,  2.18] 

rewrite  [3.7]  as  follows 

we  may 

I«(a,/3)  ■ 

•  m 

/dv/dT  bs(a,^v,T)q8(a,p,v,T) 

t3.9] 

where 

q8fa,/3,v,T)  ■ 

hs(a,»,v,T) 

[3.10] 

(v“-p-iiiv)  [(a-v)#  -  r  -  iM(a-v)] 

bs(a.p,v,T)  ■  | 

’g  «r?  .1(“-v)V<il)j//j£  S*  .‘^dxdsj 

L3.il] 

p  .  (v#+t*)^  r  .  [{a-v)*  + 


[3.12] 


HYDRONAUTICS,  Incorporated 


-15- 


and  p,  T  are  real  and  positive  for  v,T,a-v,0-T  real. 

s 

Due  to  the  use  of  [3.6  J  the  function  q  (a,j8,v,T)  and 

g 

b  !of,/3iv,T)  have  the  following  symmetry  property:  under  a 
change  of  variables  a-v  »  v,  /3-t  =  t 


qS(a,0,v,T)  =  qS(a,0,v,T) 
bS(a,^,v,T)  =  b8(a./3,v,T) 


13.13] 


The  contribution  of  qs  in  [3.9]  under  integration  over  the 
variable  v,  for  instance,  may  be  further  split  into  an  imaginary 
part  resulting  from  the  semiresidue  at  the  poles  of  q8  [3.10] 
and  a  real  part  from  a  Cauchy  principal  value.  The  poles  of 
the  quotient  l/(v“-p)  of  q  are  at  v  -  *  v  ,  where  similarly 
to  [8.22]  P 

$ 

Vp  -  [1  *  (l+*T»)*]j  t3.1*] 

and  they  have  to  be  circumvented  from  below.  The  other  two 

poles  are  associated  with  the  quotient  l/l(a-v)#-  r]  in  [3.10], 

-ut  due  to  the  symmetry  of  q*  [3.13]  it  is  enough  to  account 

for  the  semiresldue  only  at  s  v  and  to  take  the  result  twice. 

P 

Hence, 
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0» 

I§a,p)  =  Fl J  dT[bs(a,f,v  ,T)q®(a,p,v  ,T)+bs(a,£,-v  »T)Q'(a,p,-v  .t)] 


with 


cm  an 

+  /drj"  4v[bS(a,f,v,T)qs(11(p>v,T)+b5(o,d>-v,T)q8(a,d>.v,T)] 


13.15] 


q*KS,v,t> - jvh  (»rp,v,n 

p  [(a-v)*-r]  (?v*-l) 


[3.16] 


The  symbol  J-  stands  for  the  Cauchy  principal  value  at  the 
poles  of  the  Int  $rand.  Due  to  [3.133  we  have  the  following 
symmetry  properties 

}  [3.17] 

§  (ft#*** 

where  b#  is  complex  conjugate  9f  b#. 

With  the  aid  of  [3.15]  and  [3*1?]  we  may  write  the  spectral 
functions  in  their  final  form  at  fellows 
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S?0)-1C?(^)  =  211 f(a  = 


-2ir 


w 

/  lb%p,£,VjT)qJ(a  ,*,yT) 


-  b  (-a  ,0,v  >q  (-3  ,£,v  ,  ;  jdt  + 
P  P  F  p  p 


«  * 

+  2i  f  dr  JdvlbS(cio?£,v,T 


)q  (ci  }  + 

P 


+  b  '-ey0,v,nq  (-ty£sv,T|] 


C3.1B ■] 


The  function*  (3.n,  It. If,  3.18]  end  q*  [3.ie]  *r*  rfttl 
*nd  do  not  depend  ©n  the  shape  ->f  the  hip.  feg  f  3.  n  j  u  g*ner- 
complex  and  depend©  quadratic  ally  on  Cf/lx. 


3*  The  C 


b 


of  §*  and.  S«- 


■We  -consider  now  the  second  order  body  effect. .  first*.  #e 
transform  the  expression  of  I»  [iM,  i.io]  by  in  teg  retire  by 
carte  and  by  aetusing  that  along  the  contour  of  .§  f  <?  Q  f.  the 
result  being  valid  even  for  a  fiat  pet to# )  to  obtain  ' 


I#Ki 


.  *bc? 


0}  *  la^'n.r) 


*  y 


s- . 


„>Z  1®X 

e  e  sad 


1-3 


f  x,i%  2) 


“V 

/  ft*. 


$x 


m 


(x»O»0$. 


1.3*  *“ 
e  &x 


13.19 
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Hence  I'a  has  been  split  into  an  integral  over  the  center  plane 
Island  a  waterline  integral  Ia^. 

be 

We  transform  the  expression  of  la  by  substituting  for 
depx/dx  and  dcpi/dz  in  [3-193 


depi  (x,0,z) 
Sx 


2ir 


00  00 

/-/ 


Scpv 


dT  — -  e 

_os  _co 


-ivx 


[3.20] 


Scpx  (x,0,z) 


00  00 


1  f  f  .  5cpl  -ivx 
=  or  J  dv  J  dT  — :  e 

dz  -®  -»  dz 


[3.21] 


and  after  using  [3.4,  3. 5/3.1]  and  an  integration  by  parts  over 
x  we  obtain 


00  00 


j=IV 


la  G  ( a  j  )  =  /  dv  /  dT  Z  qjb 


-CO  .00 


b 


j=i 


[3.22] 


where 


y  b  av  b 

”  5  qTi  =  ttt-tt  >  9i 


_  _  YP+av  .  .  -b  _ 

a-v  '  *13  '  p(a-v)  *  ^III  "  p(a-v)  '»  4IV  " 


2(yp+av) 
(a-v)  x  v^-p-inv) 


[3.23] 
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f 

r ..  1 

*5 

-p  J  2-Z  |  =. 

-  r  sign(z-z) 

bb  =bb 

=  j+rff Jxdi^  eyzel<a-v)xf[ dial  —  elv* 
s  dx  S 

e-p|z-z| 

PP 1 2+2  1 

III  IV 

J 

l  -/ 

bl 


Next  we  consider  Ig  [3.19]  which  under  substitution  of 
[3.2l]  and  after  using  [3.1*  3.2]  and  an  integration  by  parts 
over  x  becomes 


QD  00 

Is C(aj  £)  =  y  dv  J  dT  qb 


-CO  _00 


[3.25] 


where 


b 

= 


V  (a-v)  (v“-p-inv) 


2v® 

3" 


[3.26] 


hb  -  1 

V  '  *4  7 


/  dx  ei(a“v)x]  JJ" d=d-  MMI  ePZeivx 


-t 


5x 


&X 


[3.27] 


Summarizing,  Ig  [3.19]  has  the  following  expression 
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y  co  oo 

l?(a,0)  =  £  f  dv  J  dT  q^(a, 0,v,T)bk(a, 0,v,t)  [3.28] 

0=1 

Equation  [3.28]  has  the  same  structure  as  [3.9];  we  proceed, 
therefore,  along  the  same  lines:  we  separate  the  contributions 

of  q1?  at  the  poles  v  =  [3.14]  and  on  the  rest  of  the  v  axis 

J  P  b  b 

and  use  the  symmetry  properties  of  q  and  b  which  are  similar 

J  J 

to  [3.17].  The  final  result  for  the  spectral  functions  is,  like 
in  [3.18], 

sJ(/9)-iC?(0)  =  21I?(ap,f5,v,T)  = 


J=V 


=  -4tt 


f 

Z  J  dT^(aP^>v^T)q^  - 


«j=IV 


M>PX  P‘ 


-  bq(-a  ,0,v  ,T)q.  ( -a  ,  j3,  v  ,t )  ]  + 

y  p*^  p*  'M,pv  p  P 


V  f  f  ’ 

+  41  J]  J  d tJ  dv[bj(ap,i5,v,T)q^(ap,iB,v, 


t)  + 


J=I 


o  o 


+  bj(-ap,p,v,T)q^(-ap,£,v,T)] 


[3.29] 


where 


lIV,p 


2v(yp+av) 
~(a-v) (2vm-1) 


b _ 2v3 

qV,p  “  (a-\>)  (2v*-l) 


[3.30] 
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The  principal  value  in  the  last  integral  of  [3.29]  is  at 

v  =  v  solely.  The  singularity  at  v  =  a  of  q.  [3.23,  3.26]  and 
b 

q . 


?3.30] 


is  removable  since  b .  =  O(a-v)  for  a-\r*0.  Like  in 

_ _ _  u  b  0 

[3.18]  the  function  q.  and  q  in  [3.29]  are  independent  of  the 

J  ^  <J ,  P 

shape  of  the  ship  while  depend  quadratically  upon  df/Sx. 

I  l 

4 .  The  Computation  of  Sa  and  Ca 


Last  we  consider  the  waterline  integral  second  order  effect. 
By  substituting  [3.20]  in  la  [2.2l]  and  after  using  [3«3,  3*2, 
2.18]  we  obtain 


la  (a,/3) 


co  a 

■  f»f 

-00  -CO 


dT  q^( v, t )b^(a, v, t ) 


[3.31] 


where 


l 

q  (v,t) 


2y 

v^-p-ipv 


[3.32] 


[3.33] 


By  the  same  procedure  of  splitting  the  integral  over  v  in  [3.31] 
as  in  the  preceding  sections  we  arrive  at  the  final  equation  of 
the  spectral  functions 
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SaO)-iCa(iB)  =  2ilJ(a  ,  jB)  = 


w 

=  ~hirJ  qp(^p)Cb^(ap,^vp,T)-b't(-ap>^,vp,T)]dT  + 


+  4 


1  JdrJ-dv  q4(v>T)[b*(a  ,0,vjT)+b*(-a  ,/3,v,t)] 


o  o 


[3.34] 


where 


.  2 y3 

Vv)  ~  2v8-1 


[3.35] 


The  structure  of  [3.34]  is  similar  to  that  of  [3.18]  and 
[3.29]  and  it  has  the  same  properties  as  the  other  spectral 
functions . 

5.  Summary 

Equations  [3.18],  [3.29]  and  [3.33]  render  the  various 
second  order  free  waves  spectrum  functions  in  a  compact  and 
similar  form.  The  total  spectrum  functions  are  obtained  by 
summation 


C*  ( ^)  “iSj  ( ^)  -  Cj+C| +C| -i ( S|+S|+Sj  ) 


[3.36] 
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In  the  most  general  case  in  which  the  shape  function  f(x,z) 
is  not  given  in  an  analytical  form,  two  numerical  integrations, 
over  x  and  z,  are  required  in  order  to  obtain  a  point  on  the 
curves  of  Cx  ( £)  or  St  ( /?)  [2.28,  2.18].  To  compute  the  value  of 
one  of  the  components  of  Ca  or  Sa  [3*  35 },  for  a  given  jS,  six 
numerical  integrations  are  generally  needed  namely  over  x,  z,  x, 
z,  v  and  t.  This  amounts  to  a  huge  volume  of  computations  and 
current  calculations  of  second  order  effect  for  ships  of  arbitrary 
shapes  do  not  seem,  therefore,  to  be  of  practical  interest  at  the 
present  time. 

Considerable  simplifications  are  achieved,  however,  if 
f(x,z)  has  a  sufficiently  simple  analytical  expression  or/and 
if  the  ship  is  symmetrical,,  Due  to  the  fundamental  interest  in 
computing  second  order  effect,  it  is  worthwhile  to  explore  first 
these  simple  cases. 


IV.  THE  SPECTRAL  FUNCTIONS  FOR  SEPARABLE 
SHAPE  FUNCTIONS  AND  SYMMETRICAL  SHIPS 

1 .  Separable  Shape  Functions 

If  the  function  f(x,z)  is  sufficiently  regular  it  may  be 
expressed  under  quite  general  conditions  as  a  separable  sum 

k«n 

f(x,z).  =  £  hk(x)mk(z) 

K~1 


[4,1] 
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where  trigonometric  or  polynomial  functions  have  been  generally 
used  for  representing  h^  and  m^.  Considering,  for  the  sake  of 
simplicity,  only  one  term  of  [4. l"l,  written  generically  as 
h(x)m(z),  leads  to  the  following  expressions  of  the  functions 
related  to  the  spectral  functions  of  Chapters  II  and  III: 

(i)  The  first  order  function  Ix  [2.18]  becomes 


Ii(a,0)  =  f  dx^e  f  m(z) 

J  dx  J 

-l  -t(x) 


"'ieryZ  dz 


[4.2] 


(ii)  The  functions  bS  [3.11],  b^  [3.23,  3.27,  j  =  I,  ...V) 
and  b”1,  [3*32],  which  are  the  only  ones  depending  on  the  ship 
shape,  may  be  written  as  follows 


/b8 


V, 


=  5^  / 

-l 


~  dh  _i(a- 


dx  ~  e 
dx 


a-v)x  J 


T  (a,jB,v,T,z ,z)' 

J  ar  J  <>MTj(a,/3,v,T,i,z) 

'*  ^  -t(*>  -t(S)  \t4(-.».v.t  J.i)i 


dh  ivx 

—  e 


t*.3] 


where 
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TS  = 

m(z) 

erz  m(z)  e^Z 

ii 

m(z) 

eyz  m(z)  e~^z 

~z '  signj'z-z  | 

[4.5] 

M 

II 

m(z) 

e72  m(i)  e""^ 

-I| 

[4.6] 

,b  = 

III 

Tb  : 
IV 

=  m(z)  eyZ  m(z) 

eP(z+z) 

[4-7] 

n 

* 

m(z)  m(5)  e^Z 

6(z) 

[4.8] 

and  6(z)  is  the  Dirac  function. 

Hence,  if  the  integrals  in  [4.3]  can  be  carried  out  in 
closed  analytical  form,  finding  a  point  on  one  of  the  C,(/3)  or 
S#(j3)  curves  requires  two  integrations  solely,  over  v  and  t. 
However,  due  to  the  similar  structure  of  the  b  functions,  much 
computer  time  can  be  saved  by  a  rational  programming. 

The  functions  q8  [3.10],  q*|  [3.23],  q4  [3. 32],  T8  [4.4] 
b  J 

and  Tj  [4.5  -  4.8]  are  slowly  varying  functions  of  v  and  t.  In 
contrast,  the  integration  over  x  and  x  in  [4.3]  renders  oscil¬ 
latory  functions  of  v»  the  oscillations  become  rapid  as  l  ■ 
l ' g/U ' •  becomes  small.  Accurate  numerical  integration  over  v 
of  the  integrands  of  the  spectral  functions  imposes  a  dense 
grid  in  the  latter  case,  but  considerable  computer  time  may  be 
saved  if  the  Integrals  over  t  of  the  slowly  varying  functions 
are  evaluated  by  Interpolation  from  a  sparser  grid. 
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Finally,  if  the  contour  of  S  is  rectangular,  i.e.,  |x I  <  t, 

0  >  2  >  -t  and  t  =  constant,  the  first  two  integrals  and  the 
last  two  integrals  in  [4.3]  separate  and  the  function 

L(a,v)  -£(/*!  j  j  J  A  |  6^|  [».9] 

can  be  computed  once  for  all  second  order  terms. 

2.  Symmetrical  Ships:  f(x,z)  ■  f(-x,z) 

In  this  cast  the  function  lx  [4.2]  becomes 

li(ot,0)  ■  i  ff  —  evz  sin  oix  dx  [4.10] 

and,  as  well  Known,  the  cosine  spectrum  function  C|(0)  vanishes, 
while  according  to  [2.28] 


«tO) 


sin  ax  t-Td*1 
P 


[4.11] 


As  for  the  second  order  terms,  the  functions  b8  [3.11], 
bj  [3.21,  3.27]  and  b*  [3.32]  are  purely  Imaginary.  As  a  result 
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s  s 

the  cosine  and  sine  spectrum  functions  Ct  and  SB  [3.18],  become 
equal  to  the  first  and  second  integrals  of  [3.18],  respectively, 
and  the  same  is  true  for  CB,  SB  [3.29]  and  CB,  SB  [3.34].  It 
turns  out,  therefore,  that  computing  the  cosine  spectrum  func¬ 
tions  requires  one  integration  less  than  those  necessary  for 
the  sine  functions. 

V.  APPLICATION  TO  A  SYMMETRICAL  SHIP  WITH 

PARABOLIC  WATERLINE  AND  RECTANGULAR  FRAMES 

1.  General  Formulae 

We  consider  now  the  particular  case 


f(x,z)  -  h(x)  m(z) 

h(x)  ■  l-(x#A*)  »  m's)  "1  ( | x |  <  1,  0  <  z  <  -t) 


[5.1] 


with  t  -  constant.  This  example  has  been  selected  because  It  Is 
the  simplest  conceivable  smooth  shape  and  because  detailed  mea¬ 
surements  have  been  carried  out  for  models  of  such  shape. 


For  the  first  order  free  waves  spectrum  [2.28,  4.2]  we 
obtain  immediately  the  well  known  expression 


s,(M 


-a  •  t 

8  l-«  p 

z 


sin  a l 


■  1  ■  ■  -  (cos  at- 

V 


MM  *  o 


[5-2] 
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where  a  is  given  by  [2.22]  and  Ci  by  [2.27], 

Jr 

The  spectral  functions  C*  and  [3.18]  depend  in  part  on 
the  auxiliary  functions  qs  [3.10]  and  qs  [3.16]  which  are  inde¬ 
pendent  of  the  shape.  The  function  bs  [3.11]  becomes  in  the 
present  case 

b  -  iL(a, v)  M  [5*3] 

where  L(a,v)  is  defined  by  [4,9]  and  M?  results  from  [4.3], 
Hence,  by  using  [5.1]»  we  obtain  in  the  present  case 

t(a>'))  *  -  jsfi-  tc0‘  a4  +  008  - 

■  totn  al  +  ,ln  (°-2v)0  - 

‘  >(a -vlt  t*ln  04  *  ®in  - 

'  lco*  at  '  008  fo^v)*] 

while  M*  is  given  by 
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and  are  defined  by  [3.12]. 

The  cosine  spectral  function  [3.18]  is,  therefore,  given 
by 

m 

Ct'P)  *2 rr  J  df  [L(op,vp)  MS(ap,/3,  vp#T)  qS(ap,0,\>p,T)  + 

+  Lf-Op.Vp)  M®(-ap,p,Vp,o  ts.6] 

whereas  the  sine  spectral  function  [3*18]  has  the  expression 


sf(fi)  •  -2  £  4vJ  dv  [L(op,w)  M®(op,^,v,t)  q8(ap,0,v*T)  - 

“*  o 

“  L(-ap,w)  M  t)  q  (^op>£*u»T)3  [5*73 

all  the  functions  appearing  In  [5.6]  and  [*>.7]  being  defined 
previously* 

.Similar  equations  are  obtained  for  the  body  and  waterline 
Integral  second  order  spectral  functions.  The  only  specific 
computations  are  required  by  [*.3] 
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o  o 

M^(a,0,v>T)  =  J  dz  J  dz  T^(o> ^»v>T>2jZ) 
-t  -t 


which  leads  with  the  aid  of  [4.4  -  4.8]  and  [5*1]  to 

M*  «  ±  f-  -A-  (e**>t-e’,vt)  +  4-  [l-e"^*]) 

I  p  ^  Y-p  v  ^+p  **  iJ 


[5.8] 


[5.9] 


■  ?  {?  '‘■•"'I  •  £  ^  t».») 


5.. 


My  *  M4  -  ^ 


[5.12] 


Henct,  with  bs  •  w«  obtain  th*  roliowln*  final 

expressions  for  the  spectral  functions  [3*29]  *od  [3.34] 


qJ,P(V*’VT,MJ(ap 


m 

P 


"j* 

I 

J* 

m 


L*"vV  £  q5,p("ap#l>vT)w5(Hap#IStv,,l  \ 

J*1V 


•i" 

>  > 

J  W 
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ca  • 

r 

p  - 

j=v 

S?(P)  -  dv 

{ L'Vv' 

Y  M^(a  ,  ^,m,t  )  q^(ci  ,  v,t  ) 
o  P  J  P 

O  0 

l 

LJ=I  J 

— » 

t* 
i  ii 
<  i 

_ i 

I>{  -aD»u) 

Y  ^(*ap^>v,T)q^(-ao,^>v>T) 

i 

H 

» 

*r» 

«w 

C i(fi)  -  dt  tWo  )  +  M~yvp>]M*(vp,^qp(vp}  [5.15] 

o 

S^f)  *  •fcj^d'ej  dy  tL(ap»v)  -  M-y-jtfu.e^Ns*)  [5*14} 

o  o 

2.  Discussion  of  the  Integration  Procedures 

The  second  order  spectral  functions  have  to-  be  e^puted 
nuaericaliy.  The  S*  functions  f-  .13}  and  art 

based  on  one  Integration  met  V  if  U  U  ihe  St  [$.?],  [$. It]  . 
and  [$.l6]  imply  a  double  integration.  The  »*4©r  contribution, 
in  aft  asymptotic  sense,  to  the  save  drag  [£«2&]  result?  frca*  the. 
aine  functions/  Herewith  are  a  few  observations  related  to  the  in 
teg ration  procedures  in  the  w,*  plane. 

(i)  it  is  easy  to  ascertain  that  the  behavior  at  infinity 
of  the  integrands  In  [$.6  -  5*14}  ensures  integrabUity.  the 
integrands  in  [$.4],  f$. i*]  and  [$.i$]  are  0(t'*J  for  t  -  •, 
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while  the  integrands  in  [5-7],  [5.12*]  and  [5.16]  are  0(p"3) 
for  p  =  Vv2+ts  It  is  worthwhile  to  point  out  that  if  a 

slender  body  type  additional  expansion  is  carried  out,  i .  e ., 
an  expansion  of  e  e  e  in  a  t  power  series  in  T5 .5] 
and  [5.9  -  5.12j,  the  integrals  of  the  spectrum  functions 
diverge  because  of  the  limited  validity  of  such  an  expansion 
for  large  p . 

(ii)  The  singularity  of  L(a,v)  [5-^]  at  m  =  a  is  removable. 
In  fact,  it  is  easy  to  ascertain  that  for  M~a  -  0  L(a,v)  = 

O(v-a),  which  makes  the  integrand  finite  even  for  f3.23]. 

J 

The  singularity  of  [5.9]  and  [5. 10]  at  y-p  =  0  is  also  removable. 

s 

(iii)  The  integral  in  the  \j,t  plane  rendering  S8  [5-7] 
involve  some  principal  value  computations.  As  a  matter  of  con¬ 
venience  it  was  found  that  it  is  preferable  to  carry  out  first 
integration  over  t  and  afterwards  over  v.  The  reason  is  that 
the  oscillatory  function  L(a,v)  may  be  taken  before  the  t  in¬ 
tegral,  while  the  remaining  functions  vary  slowly  with  t. 

g 

The  singularities  of  the  denominator  of  q  [3.10]  are  dis¬ 
tributed  along  a  few  curves  in  the  v,t  plane,  which  are  repre¬ 
sented  in  Fig,  2. 

Integrals  along  lines  v  --  constant  in  [5.7]  may  cross, 
therefore,  up  to  four  times  the  curves  of  Fig.  2  and  at  each  such 
point  an  evaluation  of  a  Cauchy  principal  value  is  needed.  At 
the  beginning  we  have  used  a  procedure  suggested  by  Landweber 
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(Kobus,  1967,  Appendix  1)  in  order  to  compute  the  principal 
values  contributions.  This  procedure  stipulates  selecting  the 
singularity  as  a  node  in  an  equally  spaced  grid  and  replacing 
the  integrand  by  the  derivative  of  its  regular  part  at  this 
point.  The  application  of  the  procedure  has  been  found  to  be 
tedious  in  the  present  case  for  a  few  reasons:  at  each  value 
of  v  rezoning  of  the  grid  is  necessary,  the  use  of  efficient 
numerical  integration  based  on  orthogonal  polynomial sets  pre¬ 
cluded,  the  vertices  of  the  curves  in  Fig.  2  require  special 
treatment  and  numerical  computation  of  the  derivative  is  time 
consuming. 

The  procedure  used  instead,  which  obviously  is  not  the 
only  possible  one,  was  as  follows:  in  an  integral  of  the  type 


00 


(v8-p)  [(a-v)--r] 


dT 


as  encountered  in  [5-7]>  the  singular  part  has  been  extracted 
and  I  has  been  replaced  by  I’,  where 


F(v,r)  _  _  F’(y)  F"(y)  F'"(y) 

8  p)  l(a'“'J)  S_T]  T*-T  ’  *  Tp-T'  1  *  ~  Tft-T  '  1,8 


17] 


[5.18] 
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P'(v)  =  - 


2ya  Ffv.r1) 
(a-v)a-r(a,  ) 


[5.19] 


F”(m)  =  - 


2t* 1  (a-y)8  F(v,t' 1 


[5.20] 


F'"(v)  = 


2t  1 11  (a-y)8  F(  y ,  t  * 1  1 

[Va-P(V,T"') 


[5.21] 


t*  =  v  Vva-1 


[5.22] 


t"  =  J3  +  (v-a)  V(v-a)8-l 


[5.23] 


t"  '  *  /3  -  (v-a)  V(v-a)8-l 


[5*24] 


The  integrand  in  [5.18]  is  regular  at  the  singularities  of 
Fig.  2,  including  the  vertices,  and  any  numerical  procedure  may 
be  used  in  order  to  compute  I'  provided  that  a  too  close  neighbor¬ 
hood  of  the  singularity  is  avoided. 

Similar  formulae  are  valid  for  the  part  of  s|  [5.7]  involving 

g 

q  (-a  ,  with  a  replaced  by  -a  . 

P  P  P 

(iv)  The  integrals  rendering  sf  [5.14]  and  sj  [5.16]  have 
the  only  singularities  along  t  =  t'  and  the  former  procedure  may 
be  used  in  order  to  make  the  integral  regular,  only  the  F»(v) 
[5.19]  type  of  function  being  needed. 
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VI.  PRELIMINARY  NUMERICAL  RESULTS 
1.  Numerical  Procedure 

We  present  a  few  preliminary  results  of  numerical  computa¬ 
tions  of  second  order  spectrum  functions.  These  first  results 
are  presented  on  a  tentative  basis;  they  have  to  be  furtherly 
checked  in  terms  of  accuracy  of  numerical  procedures  and  to  be 
extended  to  a  larger  range  of  Froude  numbers.  The  computer 
program  has  to  be  improved  towards  saving  time. 

The  computer  program  has  been  written  in  order  to  determine 
the  values  of  S8  [5.8],  S8  [5.14]  and  S8  [5.16].  The  computation 
of  the  C8  functions,  which  is  much  easier,  will  be  carried  out 
in  the  future. 

As  stated  previously,  numerical  Integration  has  been 
carried  out  first  along  lines  v  =  constant  in  the  v,t  plane. 

A  48  points  Legendre  procedure  has  been  applied  repetitively  to 
intervals  of  At  =  0.25.  The  integration  has  been  stopped  if  the 
increment  of  the  integral  on  such  an  interval  was  lesser  than 
10"4  of  the  sum  on  previous  intervals  or  after  25  intervals  if 
the  criterion  was  not  met.  In  the  latter  case  integration 
continued  over  25  enlarged  new  intervals  until  the  convergence 
criterion  was  satisfied.  When  the  integration  was  stopped,  the 
rest  of  the  integral  on  the  remaining  t  interval  up  to  t  =  • 
was  evaluated  by  a  24  points  Laguerre  procedure.  The  change  of 
At  from  0.25  to  0,5  showed  very  little  change  in  the  values  of 
the  integrals. 


HYDRONAUTICS,  Incorporated 


The  integration  over  v  has  been  carried  out  with  a 
Simpson  procedure,  the  v  interval  being  w/80l ,  selected  such 
that  enough  points  should  be  contained  in  a  "wavelength"  of 
the  trigonometric  functions  which  are  present  in  L(a-v)  [5.4]. 
The  integration  has  been  stopped  when  the  absolute  value  of 
the  increment  was  smaller  than  lO-3  of  the  previous  sum. 

Summarizing,  the  program  admits  l  =  t'g/U'*  and  t  =  t'g/U 

g  “k  J 

as  input  parameters  and  computes  SB,  SB  and  Sa  for  any  given 
value  of  0. 

2.  Numerical  Results  for  Given  Draft  and  Froude  Numbers 

The  first  runs  have  been  carried  out  for  t'/l'  =  0.3  and 

for  i  =  i ' g/u ' 3  =5  (F  =  0.316)  and  t  =  12.5  (F  =  0.200). 

**  sn  v  i 

The  values  of  the  sine  spectrum  functions  Sj ,  SB,  SB  and  SB 
have  been  computed  for  0  **  0,  0.5,  1,  1.5,  2,  3>  4,  5,  6.  In 
Figs.  3  and  4  we  have  represented  the  functions 


oS,b,<t 

b» 


*»o  »SBS,bj/' 
_ E _ - _ 

V"1 


[6 


which  permit  to  rewrite  the  equation  of  the  free  waves  profile 
[2.27]  in  the  more  convenient  form 
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where  e  =  b ' / l'  is  the  usual  slenderness  parameter. 

Although  conclusions  may  be  drawn  only  on  a  provisory  basis, 
it  is  worthwhile  to  discuss  a  few  qualitative  features  of  the 
curves  of  Figs.  3  and  4: 

(i)  For  F^  =  0.316  (Fig.  3a)  the  free  surface  correction, 
the  body  correction  and  the  waterline  integral  correction  are  of 
the  same  order  of  magnitude.  For  the  smaller  Froude  number 

F^  =  0.2  (Fig.  4a)  the  free-surface  effect  seems  to  tend  to 
dominate  the  other  two  contributions. 

(ii)  For  Fn  =  0.316  (Fig.  3b)  S9  =  SB  +  sj  +  Sg  is  roughly 
equal  and  of  opposite  sign  to  .  Hence,  the  second  order  term 
diminishes  the  amplitude  of  the  first  order  free  waves  and  of 
the  corresponding  wave  resistance.  For  F^  =  0.200  (Fig.  4b)  SB 
is  again  of  opposite  sign  to  Sv,  but  the  amplitude  at  the  first 
peak  of  the  Sg  curve  is  roughly  four  times  larger  then  that  of 

£>1 . 

(iii)  The  first  order  solution  is  valid  only  if  7  «  1, 
let  say  7  <  0.1  for  Fn  =  0.316  and  7  <  0.025  for  Fn  =  0.2,  i.e., 
as  well  known  not  only  7  but  c  »  7<.  has  to  be  small  as  second 
order  effects  become  more  important  at  low  Froude  numbers. 

VII.  SUMMARY  AND  CONCLUSIONS 

A  systematic  procedure  of  computing  the  second  order 
spectral  functions  of  the  free  waves  generated  by  a  thin  ship 
has  been  derived.  The  procedure  has  been  applied  to  the  case  of 
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a  ship  with  parabolical  waterline  and  rectangular  frames.  Pre¬ 
liminary  numerical  results  indicate  that  the  free  surface,  body 
and  waterline  integral  second  order  corrections  are  of  the  same 
order  of  magnitude  and  neglection  of  part  of  them  in  order  to 
derive  inconsistent,  but  simpler,  solutions  does  not  seem  to  be 
justified.  This  conclusion  has  to  be  strengthened  by  extending 
the  computations  to  various  shapes  and  a  large  range  of  Froude 
numbers. 

The  accuracy  of  the  numerical  program  is  checked  at  present. 
Later,  cosine  and  sine  spectral  functions,  as  well  the  associa¬ 
ted  wave  resistance,  will  be  computed  for  various  Froude  numbers, 
such  that  drawing  more  definite  conclusions  will  become  possible. 
These  future  results  will  be  reported  elsewhere. 
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«OUtt  3  -  THE  SINE  SFECTRUM  FUNCTIONS  OF  THE  FREE  WAVES  FOR 

F  •  0.316  (t*  5) 
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FIGURE  4  •  THE  SINE  SPECTRUM  FUNCTIONS  OF  THE  FREE  WAVES  FOR 

F  *  0.200  ( t  *  12.5) 
n 


